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Abstract It was found that the advection-dominated accretion 
flow (ADAF)+thin disk model calculations can reproduce the 
observed spectral energy distributions (SED) of the two low- 
luminosity active galactic nuclei (AGNs), provided they are ac- 
creting at ^ 0.01 — 0.03 Eddington rates and the thin disks are 
truncated to ADAFs at ~ 100 i?s {Rs is the Schw arzschild ra- 
dius) for M81 and NGC4579 dOuataert et all Il999h . However, 
the black hole masses adopted in their work are about one order of 
magnitude lower than recent measurements on these two sources. 
Adopting the well estimated black hole masses, our ADAF+thin 
disk model calculations can reproduce the observed SEDs of these 
two low-luminosity active galactic nuclei, if the black hole is ac- 
creting at 2.5 X 10^"* Eddington rates with the thin disk truncated 
at Rtr = 120i?s for M81 (m = 3.3 x lO'^ and Rtr = 80i?s 
are required for NGC4579). The transition zones with tempera- 
ture from the thin disk with 10^"^ to 10^~i° K in the ADAF 
will inevitably emit thermal X-ray lines, which provides a use- 
ful diagnosis on their physical properties. The observed widths of 
the thermal X-ray iron lines at ~ 6.8 keV are consistent with the 
Doppler broadening by the Keplerian motion of the gases in the 
transition zones at ^ 100i?s- We use the structure of the transition 
zone between the ADAF and the thin disk derived by assuming 
the turbulent diffusive heat mechanism to calculate their thermal 
X-ray line emission with standard software package Astrophysical 
Plasma Emission Code (APEC). Comparing them with the equiv- 
alent widths of the observed thermal X-ray iron lines in these two 
sources, we find that the turbulent diffusive heat mechanism seems 
to be unable to reproduce the observed thermal X-ray line emis- 
sion. The test on the evaporation model for the accretion mode 
transition with the observed thermal X-ray line emission is briefly 
discussed. 
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1 INTRODUCTION 

It is well accepted that many astrophysical systems are powered by black hole 
accretion. The standard thin disks (i.e., geometrically thin and optically thick ac- 
cretion disks) can succ essfully explain most observ ational features of the black 
hole accretion systems dShakura & Syunyaev[ll973h . The ultraviolet/optical con- 
tinuum emission observed in luminous quasars is usually attributed to the ther- 
mal radiation fr om the standard disks (SD) surrounding the massive black holes 
in quasars (e.g.. lSun & Malkanill989h . However, the standard disk model is un- 
able to reproduce the spectral energy distributions (SED) of many sources (e.g., 
Sgr A*) accreting at very low rates, and the advection dom inated accretion flows 
(ADA F) were suggested to be present in these sources dNaravan&YlL 11994 
fl995l) . In the ADAF model, most released gravitational energy of the gases in 
the accretion flow is converted to the internal energy of the gas, and the ADAF is 
hot, geometrically thick, and optically thin. Only a small fraction of the released 
energy in ADAFs is radiated away, and their radiation efficiency is th e refore 
significantly lower than that for standard thin disks (see iNarayan et all Il998l 
for a review and references therein). The ADAF model can successfully repro- 
duce the observed SEDs in many black hole sy stems accreting at low rates (e.g., 
iLasota et al.lll996HNaravan et al.lll99d.ll998l) . 

There is a critical accretion rate merit, above which the ADAF is suppressed 
and a standard thin disk is present. The ADAF can co-exist with the standard thin 
disk, when it is accreting at rates slightly lower than the critical value rhcrit- In 
this case, the ADAF is present in the inner region near the blac k hole and con- 
nects to a standard thin disk at a certain tran sition radius J?tr (see Naravan et all 



Il998l for a review and references therein). iNaravan et al.l ( Il996 ) proposed the 



ADAF+standard thin disk systems to model the obser ved spectr a of th e black 
hole accretion systems with moderate accretion rates. iNaravanI (Il996h found 
that many different spectral states observed in black hole X-ray binaries can 
be understood as a sequence of ADAF+thin disk models with varying m and 
Tti (where rti is the outer radius of ADAF in units of Schwarzschild radius, 
i?s = 2G'Mbh/c^ )- This scenario was explored in more detail by several authors 
dEsin et al.ll 19971 [1998), and they found that the different spectral states of X- 
ray binaries (i.e., the quiescent state, low state, intermediate state or high state) 
correspond to different accretion rates. It is believed that standard thin disks (or 
slim disks) are present in luminous active galactic nuclei (AGN), while ADAFs 
are present in those low-luminosity AGNs accreting at relatively low rates. The 
ADAF+thin disk systems are required for mo deling on a variety of observations 
of AGNs accreting at the moderate rates (e.g.. lOuataert et all [19991 ICaoll2003[) . 

The physical mechanisms of the transition from a thin disk to an ADAF and 
its structure are still quite uncertain, though a few different scenarios were sug - 



gested (e.g.,lMever & Mever-Hofnieisteiill994l;lHonmai[T996HLiu et al. l [T99 



iiviever & Mever-Hofmeiste" ( 1994t) proposed an "evaporation" mechanism, in 



which the disk is evaporated as heated by electron conduction from a hot corona, 
and then a quasi-spherical hot ac cretion flow is formed. An alternative scenario 
was suggested bv iHonnial (Il996l) . in which a "turbulent diffusive heat" mech- 
anism is employed to explore the structure of the ADAF+thin disk system. 
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This model was further developed, an d the global structure of the ADAF+thin 
disk system was derived numerically (" Manmoto & Katd l2000t iManmoto et all 
UOOO). All of these model calculations show that a rapid decrease of temperature 
occurs within a narrow zone between the inner ADAF and the outer standard thin 
dis k at ~ rtr- 

IOuataertetalJ(ll999l) showed that the optical/UV to X-ray emission detected 
from the nuclei of M8 1 and NGC 4579 can be well explained by an optically 
thick, geometrically thin accretion disk extends down to 100 Rs (inside of 
which an ADAF is present), provided their accretion rates are to ~ 0.01. The 
optical/UV bumps observed in these two sources can be attributed to the thermal 
emission from the outer standard thin disks, while their hard X-ray emissions 
are dominantly from the inner ADAFs. In their model calculations, the same 
central black hole mass Afbh = 4 x 10^ Mq is adopted for these two sources, 
which are about an order o f magnitude underestimated compared with the recent 
estimates dPevereux et all 1200 3: Barth et al., 200 ij). The iron Kg l ines w ere ob- 
served in these two sources lTshisaki et al.i,il996ciTerashima et al.lll998l) . In ad- 
dition to narrow iron Ka lines at 6.4 KeV, the broad line co mponents centered 
at Eko = 6.79 KeV were also observed in these two sources jPewangan et all 
12004) . The therm al X-ray line emission i s a useful diagnosis on black hole accre- 
tion sy stems(e.g., rNaravan & Ravmondill999l : rXu et am2006l) . lDewangan et al.l 
(|2004|) found that the observed lines from these two sources are too broad to be 
the thermal line emission from the host galaxy, and they suggested that these 
broad emission lines may probably be from the transition zones between the 
ADAFs and the outer standard thin disks. The temperature of the inner edge of 
the thin disk is around 10'* K, and therefore the temperature of the transition zone 
extends from ^ 10^ to ^ 10^~^ K while connecting to the ADAF. Thus, the ther- 
mal X-ray line emission can naturally originate from such transition zones. The 
observed widths of the broad lines are roughly consistent with the Dopplerbroad- 
ening by the Ke plerian motion of the gases in the transition zones at 100 Rs 
dOuataertet al.Lll999,) . 

In this work, we re-investigate the ADAF+thin disk systems in these two 
low-luminosity active gal actic nuclei, NGC 457 9 and M81, a dopting the black 
hole masses estimated by lDevereux et all (l2003h : lBarth et all f2001.) . The ther- 
mal iron Ka line emissions from the transition zones of the ADAF to the thin 
disk in these two low-luminosity active galactic nuclei, NGC 4579 and M81, are 
calculated with the transition mode l based on the turbulent diffusive heat mech- 
anism developed bv lHonmal ( Il996h . 

2 THE ADAF+THIN DISK MODEL 

The observed SEDs of two low-luminosity AGNs, MS 1 and N GC4579, were re- 
produced with the AD AF+ standard thin disk model by Ouataert "et al.l (Il999l) . In 
this ADAF+standard thin disk model, no smooth physical connection between 
the ADAF and the outer standard thin disk at rn has been included. The re- 
sulted spectra are si mply the c ombination of the spectra from the ADAF and 
the outer thin disk. Ouat aert et alj (Il999 )'s model calculations showed that the 
best fits to the observed SEDs of these two sources requke the transition radius 
Ttr — 100 for both sources. The accretion rates are to = 0.03 and m = 0.01 
for NGC4579 and M81, respectively. The viscosity parameter a — 0.1, the ra- 
tio of gas to total pressure /3 = 0.9, the fraction of the released energy directly 
heating the electrons S — 0.01, are adopted in their calculations. In their model 
calculations, the same central black hole mass Mbh = 4x10^ Mq is adopted 
for these two sources, which are about an order of magnitude underestimated. 
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The masses of these two black holes were est imated as Mbh — 7 x 10^ Mp, for 
M81, and Afbh 5 x 10^ Mq for NGC4579 (lDevereux et all 120031: iBarfli et all 
[2001). Thus, we have to re-calculate the global structures of the ADAFs sur- 
rounding t he black holes i n these two sources. We employ the approach sug- 
gested by Manmotol (l2000 l) to calculate the global structure of an accretion flow 
surrounding a Schwarzschild black hole in general relativistic frame. All the ra- 
diation proce sses are in cluded in the calculations of the global accretion flow 
structure (see iManmotol [2000. for the details and the references therein). It was 
pointed out that a significant fraction of the viscously dissipated energy could 
go into electrons by magnetic re connection, if the magnetic fields in the flow are 
strong dBisnovatyi-Kogan & Lovelac e. 1997, 2000). They argued that 5 can be 
as high as 0.5 . Ther efore, we adopt a conventional value of 5 — 0.3 (e.g., 
IWu. Yuan&CaoLl2007l) . and tuning the accretion rates m to fit their observed 
optical/UV/X-ray continuum spectra. The irradiation of the outer thin disk by 
the ADAF is included in our calculations. We include an empirical color correc- 
tion given bv' Chiand (l2002h for the thermal emission from the outer thin disk, 
which can reproduce the non-LTE disk spectral model of iHubenv et al.l (1200 lb . 
The flux due to viscous dissipation in the outer region of the disk is 

which is a good approximation for i?ti ^ Rm- The outer thin disk may be 
irradiated by the incident photons from the inner ADAF region. The flux due to 
the irradiation in the outer thin disk Fi^^{R) can be calculated by assuming the 
incident photons are reprocessed as thermal radiation f rom the thin disk, while 
the global structure of the inner ADAF is available (see ICao& Wang. .2006. for 
the details). The local disk temperature of the thin cold disk is 

rj. .p. [F,,,{R) + F,„{R)Y" 

Jdisk(-r£j = , (^) 

by assuming local blackbody emission. In order to calculate the disk spectrum, 
we include an empirical color correction for the disk thermal emission as a func- 
tion of radius. The correction has the form (Chiang, 2002) 

f irj. ._r_ (/0O~l)[l + eXp(-l^/A^)] 

/coH-'diskJ Joo , r/ \ / A 1 ' ^ J 

1 + exp[(i^p - v\y)//\v\ 

where Vp = 2.82kBTdis]i/h is the peak frequency of a blackbody with temper- 
ature Tdisk- This expression for /coi goes from unity at low t e mpera tures to /oo 
at high temperatures with a transition at i/b ~ t'p- IChiand (l2002l) found that 
/oo = 2.3 and vt, = = 5 x 10^^ Hz do a reasonable job of reproducing 
the model disk spectra of lHubenv et al.l (1200 lb . The disk spectra can therefore be 
calculated by 

„ oo 

^ I c2 j c2 7 /,4,i[exp(W/coifcBTdisk)-l]' ^' 



'■tr 



We find that the best fits to the SEDs of these sources require: to = 3 x 10 ^ 
and rtr = 120 for M81; to = 3.3 x lO'^ and n,. = 80 for NGC4579. We find 
that the observations in optical/UV and X-ray bands of these two sources can be 
fitted quite well (see Fig.[T]i- 
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3 TRANSITION ZONE OF THE THIN DISK TO THE ADAF 

The global structure of the ADAF+standard thin disk systems based on the as- 
su mption of addi ti onal turbulence visco sity in the transition zone was derived 
bv lHonm"al (Il996l) . iManmoto et alj (I200Q.) performed numerical calculations on 



the ADAF+thin disk system based on lHonmal(ll996h 's model, and obtained the 



global structure of the ADAF smoothly connecting to the outer thin disk. They 
extended the Honma's analytical solution to be able to describe the global struc- 
ture of the ADAF+thin disk system in cluding the transition zone quite well com- 
pared with the numerical results (see iManmoto et all 1200 0. for more details). 
The Honma's analytical solution for the ADAF+standard thin disk system is 
only valid for fully advection dominated flows (/ = 1), where / is the ratio of 
the advected energy to the viscously dissipated energy in the ADAF. By intro- 
ducing a new paramet er a, the extended Honma's analytical solution derived by 
IManmoto et alJ (l2000l) can deal with the ADAFs with partial cooling (/ < 1), 

V = - {R/Rt,r\vK{R). 



n = 



^m^iR/R^^r/^n^iR), 



^[l~iR/Rur]'^'vKiR). (5) 
The relation between / and a is 

(a + 5)(3-a)2 ' ^' 

As we are focusing on the transition zone, ax = is adopted, because the struc- 
ture of the transition zone is almost independent of the value of this par ameter 
(corresponding to the turbulent energy transport) (lManmoto et all l2000l) . Their 
derived global structure (either the analytical one or that derived numerically) is 
one temperature (i.e., the temperature of electrons is the same as that of ions), 
which prevents us from calculating the spectrum of the accretion flow directly. In 
this work, we include all the radiation processes in our calculations for t he global 
structure of the inner ADAF, while we only use IManmoto et all (l2000l) 's model 
to connect our ADAF solution to the outer thin disk. We assume the efficiency 
of the released gravitational energy 77 = 0.1, and a fraction / of it is advected in 
the flow. We have 

^ADAF = - /)A^EddMbhC^, (7) 

which can be derived by integrating their observed continuum emission (sub- 
tracting the optical/UV bumps from the outer thin disk regions). Thus, the value 
of / can be estimated if the accretion rate m is known. The value of / can also be 
derived from our calculation on the global structure of the inner ADAF, because 
the radiation processes have been considered in the global solution of the ADAF. 
We find that the value of / estimated with Eq. O is consistent with that derived 
from the global solution of the ADAF. In order to explore how the thermal line 
emission varies with the accretion rate m, we loose the constraints on the accre- 
tion rates to, i.e., allowing m to vary about an order of magnitude from the values 
derived from the model fits on the SEDs of these two sources, in our calculation 
of the thermal line emission from the transition zone. We can estimate / as a 
function of the accretion rate to, provided Ladaf is derived from the observed 
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SEDs. The value of a can be calculated with Eq. ^ after / is known. The struc- 
ture, and then the thermal line emission, of the transition zone can be calculated 
as a function of the accretion rate m. Using Eq. Q and (|6]l, the structure of the 
transition zone is available for calculating its thermal X-ray line emission as a 
function of accretion rate. 

4 THERMAL X-RAY LINE EMISSION FROM THE TRANSITION 
ZONES 

Besides the narrow iron Ka lines at 6.4 KeV, the broad line components centered 
at EKa = 6.79 KeV with equivalent width EW=287 eV for NGC4579, and 
EW=1 01 eV for M81, were observed with XMM-Newton by iDewangan et all 
(12004'). The widths of these two lines with Gaussian fittings are a — 231eV and 
188 eV for NGC4579 and M81, respectively. 

The thermal X-ray line emission can be calculated, when the physical prop- 
erties of the plasma, i.e., the temperature, density and metallicity are specified. 
Using the models of the ADAF+thin accretion disk systems described in the 
previous section, we can calculate their thermal X-ray line emission from the 
accretion flows surrounding the black holes in these two low-luminosity active 
galactic nuclei, NGC 4579 and M8 1 . In the ADAF+standard thin disk system, 
the inner ADAF is very hot (Te ^ 10^~^''K) and the outer thin disk is cold 
(Te ~ 10'*~^K). Thus, the inner ADAF is too hot, the plasma is almost com- 
pletely ionized, to produce thermal X-ray line emission, while the outer thin disk 
is too cold for thermal X-ray line emission. For our interested H-like and He-like 
iron line emission, they are most probably emitted from the transition zone of 
the ADAF to the outer thin disk, because the temperature of the transition zone 
is between - lO^-^ - IQS-io K. 

The line luminosity iuno can be calculated by integrating over the transition 
zone, 

Lunc=/" ne{rfeunc{r)^—^-^^iTTrH{r)Rldr, (8) 

Jrtr — Sr ''"line!,''] 

where rie is the electron number density, Sr is the width of the transition zone 
connecting the inner ADAF to the thin disk, rn- is the inner radius of the thin 
disk, H — Cs /il is the vertical half-thickness of the transition zone, and 
TUnc ~ 2iJfcii„e is the optical depth of the emission line in the vertical direction 
(feline is the line absorption coefficient). The line emissivity eiino('') as a function 
of temperature is calculated wit h the standard softw are package Astrophysical 
Plasma Emission Code (APEO dSmith et all l200lh . The APEC code includes 
collisional excitation, recombination to excited levels and dielectronic satel- 
lite lines(see Smith et al .. 2001, for the details). It ignores photo-ionization, 
which is not importan t in these two low-luminosity AGNs (see discussion in 
iDewangan et al.ll200 4'). We assume ionization equilibrium in the plasma. In the 
transition zone, the density /temperature of the the flow changes by orders of mag- 
nitude from the inner optically-thin ADAF to the outer optically-thick thin disk. 
The optical depth of the flow can be very large in the outer part of the transition 
region, and the absorption and radiative transfer of the line is considered in our 
calculation. The X-ray continuum emission of the ADAF+thin disk system is 
dominated by the bremsstrahlung and the Comptonization of the soft photons in 
the ADAF. Considering energy resolution of the observations on the X-ray line 
emission is quite low, we only calculate the total equivalent widths of H-like and 
He-Uke iron lines. 
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Fig. 1 The ADAF+standard thin disk spectral models for M8 1 and 
NGC4579 (solid lines). The dotted lines represent the spectra of ADAFs, 
while the dashed lines are the spectr a of the out er standard thin disks. The 
observed data points are taken from IHoI(I1999I) . 



In Fig. 121 we plot the thermal X-ray line luminosity as functions of accretion 
rate ifi based on the transition model given by Manmoto et al. (2000). Besides 
the solar metallicity, we also calculate the X-ray line emission from these two 
sources for five times solar metallicity for comparison. 

5 DISCUSSION 

Our best fits to the observed continuum spectra of these two sources show that 
the outer thin discs are truncated at rtr = 120 and 80, for M81 and NGC4579, 
respectively. The two important disk parameters, m and rtr are mainly deter- 
mined from the comparison of the outer disk spectra with the observed opti- 
cal/UV continuum emission in these two sources, and the observed X-ray con- 
tinuum spectra can then be naturally reproduced by the ADAF models. The ob- 
served widths of the broad lines are roughly consistent with the Doppler broad- 
ening by the Keplerian motion of the gas in the transition zone at 50-150 i?g 
dbewangan et al.>j2004) . The ratio of the transition radii of these two sources is 
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Fig. 2 The total H-like and He-like iro n line luminos ities as func- 
tions of accretion rate m calculated for the lHonmal (Il996b 's model. The 
solid lines correspond to the cases with the solar metallicity, while the 
dashed lines are for the cases with five times solar metallicity. The ob- 
served values for these two sources are marked in the figure (dash-dotted 
lines). The dotted lines represent the accretion rates derived from the 
ADAFh-SD model fitting on the observed SEDs. 



1.5, which is consistent with the observed line width ratio ~ 0.814, because the 
width a oc r'^^^'^ ■ 

Comparison of ou r calculations based on the transition zones given by 
iManmoto et al.l ( |200C /) (see Fig. |2] and §3 and 4 for the detailed description) with 
the observations show that the accretion rates m ^ 7.1 x 10^^ and 6.5 x 10^* 
for M81 and NGC4579, respectively (see Fig.|2]l. These seem to be inconsistent 
with the optical/UV spectra observed in these two sources, which requires the 
accretion rates to be m = 2.5 x 10"'^ (M81) and 3.3 x 10"^ (NGC4579), unless 
the transition radii deviate significantly from ^ 100 (e.g., an order of magni- 
tude smaller than 100). However, the widths of the thermal X-ray lines provide 
strict constraints on the transition radii, if the thermal X-ray line emission does 
originate from the transition zones. The dependence of our results on the metal- 
licity can be understood with Eq. dHJ. A larger metallicity will increase the line 
emissivity and line optical depth at the same time, while the line luminosity and 
EW will slightly increase with a factor 1 — e^^"""'. Of course, the structure of 
ADAF+thin disk systems based on this model is rather simplified, in which only 
one temperature is considered (i.e., the electrons have the same temperature as 
the ions). Thus, we cannot rule out this model only from their thermal X-ray line 
emission. In the transition zone, the density is very high so that one temperature 
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flow is a good approximation, as the Coulomb interaction between ions and elec- 
trons in this region is very efficient. Our calculations of the thermal X-ray line 
emission from the transition zone have not been affected by this one temperature 
assumption. The more detailed calculations on the ADAF+thin disk systems in- 
cluding energy equilibrium between electrons and ions in the accretion flows 
may help to test this model, which is beyond the scope of this work. 

An alternative model for the transition o f a thin disk to an ADAF is the "evap - 
oration" mechanism initially suggested by Meyer & Meyer-HofmeisteJ d 19941) . 
in which the thin disk is evaporated as heated by electron conduction from a hot 
corona, and then a quasi-spherical hot accretion flow is formed. There is a very 
thin layer between the cold thin disk and the hot corona, of which th e electron 
temperature may v ary from ^ 10®^ K to '--^ 10^ K in the corona (e.g. lLiu et al.L 
Il995lll997ll2002l) . Such a layer, of course, will emit thermal X-ray line emis- 
sion. In principle, the evaporation model for the transition of a thin disk to an 
ADAF can be tested by the observed thermal X-ray line emission. However, 
the structure of this disk corona system is complicated and is only av ailable by 
numerically integrating a set of ordinary differential equations (e.g., iLiu et al.l 
l2002h . which prevents us from testing this model in this paper. 
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